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1. Introduction 



Two Higgs doublets models (2HDM) is one of the simplest, non-trivial extensions of the 
standard model Higgs sector, which appears for example in the Minimal SuperSymmetric Model 
(MSSM). For a general introduction to 2HDMs and the Higgs sector of MSSM we refer to [jl|, 0. 
A general 2HDM can therefore serve as an effective theory when searching for physics beyond the 
standard model. In this contribution we describe the program 2HDMC [g] and how it can be used 
to explore the physics of general CP-conserving 2HDMs 1 . 



2. General two Higgs doublet models 



2.1 General two Higgs doublet model potential 

We consider the potential for a general model with two complex SU (2)l doublets with hyper- 
charge Y=l: <I>i ,<t>2 which is invariant under global SU (2) transformations, <I> a — > Uab'&b, and that 
is gauge invariant and renormalizable 



m^ 2 <t>\<t> 2 + h.c. 



2" 1 V" 1 " V ' 2 
-f<{iA 5 (<I>l<I> 2 ) 2 + 



Ae )+h [&l<h 



<t>\<t>2 I +h.c. 



Demanding that the potential is real implies that m\ Y , m\ 2 , Ai_4 are real whereas m\ 2 , ^5-1 in 
general can be complex. However, if we demand that there is no explicit CP-violation then the 
latter parameters also have to be real. 

2.2 Electroweak symmetry breaking 

As usual the electroweak symmetry is broken by non-zero vacuum expectation values (vev) 
of <l>i and/or <I> 2 . After applying the minimization conditions of the potential the m 2 n and m\ 2 
parameters can be traded for the vev's of the two doublets: vi = vcos/3 and v 2 = ve'^ sin/3 with 



(V2G F 



-1/2 



246 GeV. Here ^ is a possible phase that we put to zero by demanding that 



there is no spontaneous CP- violation. The two doublets can then be written as 



1 / A/2(G + cosj3-// + sinj8) 
a/2 I vcosjS — /jsina + //cosa + i (G°cosj8 — Asinj3) 



<1>2 = — T= 



V2(G + sin/3+// + cos/3) 



a/2 I vsin/3 +/zcosa + //sina + i (G sinj8 +Acosj3) / 

where G + and G° are the Nambu-Goldstone bosons that give masses to the W and Z bosons respec- 
tively, H + is the charged Higgs boson, A is a CP-odd and h, H are CP-even neutral Higgs bosons 
(with nih < mu), and tan /3 defines a basis in <I> space. (Also note that there is a special basis, the 



'Admittedly, given the format for these proceedings, this description will have to be quite limited and incomplete 
and we apologize for that. For a more complete description we refer to p|| and references therein. 
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Higgs basis, where only one of the fields develops a vev and therefore tan /3 = or cot /3 = 0). With 
these conventions the couplings of Higgs and electroweak gauge bosons are given by the invariant 
s B-a = sin(/3 — a) and cp_ a = cos(/3 — a). Finally we note that in addition to the above parame- 
terisation of the potential it can also be useful to use the masses of the Higgs bosons as parameters. 
Thus in the 2HDMC program the user can choose between the so called general parameterisation, 
{m\ 2 , A1-7, tan/3} or the physical one {m\ 2 , tnh, mu, «U, m H ±, sp_ a , Xe-i, tan/3} in addition to 
the Higgs basis which is the same as the general one but with tan /3 = and m H ± replacing m 2 l2 . For 
backward compatibility the program can also handle the basis used in the Higgs Hunter's Guide [jj]]. 
It is important to note that the program uses conventions where tan/3 > and sp_ a > 0. 

2.3 Possible additional symmetries 

In addition to the symmetries already discussed it is also possible to impose other interesting 
symmetries on the potential. First we note that one can demand an additional U{\)pq (Peccei- 
Quinn) symmetry [Q] by setting m\ 2 = and As_7 = 0. The discrete version of this symmetry is 
to demand an exact Z2 symmetry such that the potential is symmetric under $1 — >• 3>i, <&2 — > —^2 
which leads to m\ 2 = and A6-7 = 0. This symmetry is softly broken if m\ 2 / 0, which leads to 
flavour changing neutral currents beyond tree-level. Finally in the case of a supersymmetric theory 
the potential takes a special form which at tree-level leads to 

g 2 + g a g 2 -g' 2 g 2 2 2 

X\ = A2 = — , A3 = — , A4 = — — , ^5=^6=^7 = 0, m l2 = m A cos /3 sin /3 . 

For convenience the program contains two special methods for setting up: a tree-level MSSM 
model with parameters and tan/3, as well as a so called inert doublet (ID) model [||] with {m| M , 
m 1 ®, ra^ D , ra^ ± , A2, A3 } as parameters. 



3. Theoretical constraints 

The 2HDMC program also contains methods for checking theoretical constraints on the pa- 
rameters describing the potential, namely positiviy, perturbativity and tree-level unitarity. 

Demanding that the potential is bounded from below one has the constraints [f7|] : X\ > 0, 
A2 > 0, A3 > — VA1A2. In addition, in the case X(, = Xj = one has A3 + A4 — |As| > — VA1A2 
whereas if X&, A7 7^ one gets A3 + A4 — A5 > —\JX\Xi as well as more complicated constraints 
which take too much room to include here, but are included in the program. 

The constraints from perturbativity are obtained by considering the cross-sections for 2 — > 2 
Higgs scattering processes, which can be written as a perturbative series in the square of the corre- 
sponding quartic couplings X^ jHjHkHi /(l6n 2 ). In order for this series to make sense the couplings 
cannot be too large and the program contains a method for checking this with the default limit 
being X HiHjHkHl < 4n. 

A similar constraint on the parameters of the Higgs potential is obtained by requiring tree-level 
unitarity [||] for HH and HVl scattering (with Vl denoting the longitudinal parts of W or Z), which 
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corresponds to putting limits on the eigenvalues Aryp) of the following scattering matrices 



\6kS, 



(2,1) 



167T5 1 , 
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The default value for the limit is |A, 



< 16k. 



4. Yukawa sector 

The program contains a general Yukawa sector parameterized in terms of the matrices p L , p D , 
and p u in flavour space for leptons, down-type, and up-type quarks respectively: 

-J2V =D P ' P Dh + D p J. p DH + iDy 5 ^-=DA 

V2 V2 V2 

pay p a uh v p_a_H__p_a UH _ {JJ F m 



+ 



U{V ckm P D Pr - p u V CKM P L }DH + + vp L P R LH + + h.c. 



where P R / L = (1 ± 75) /2 and fc f = \flM F jv are diagonal matrices determined by the fermion 
masses. The only constraint on p F is that they are symmetric. 

Models with non-diagonal p F will lead to non-minimal flavour violating charged currents as 
wells as flavour changing neutral currents. This can be avoided [^] by imposing a Z2 symmetry on 
<J>i, <J>2 and Ur, Dr, Lr such that each fermion type only couples to one Higgs doublet. In turn this 
leads to the restriction p F = K F cot/3 or p F = — K" F tan/3, which means that there are four different 
types of Z2 symmetric 2HDMs. For simplicity the program also contains methods for specifying 
the Yukawa sector by simply choosing the type. 

5. Partial decay widths 



With both the Higgs sector and the Yukawa sector being specified the program can then be 
used to calculate partial decay widths and branching ratios for range of processes including: 
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• H-tff with optional (N)LO QCD corrections 

• H — ^ gg with optional LO QCD corrections 

• H-^HH 

• H — > HV* including off-shell vector bosons 

• H — > VV* including off-shell vector bosons 

• //—)• 77 

• t->H + b 

where H = {h,H ^A,^} denotes any applicable Higgs boson and V = {Z, W}. 



6. Experimental constraints 

Before using the program to make predictions in a specific model one may also want to take 
into account existing experimental constraints. This is not a simple task and there are dedicated 
programs for this. Therefore, the program has been limited to having methods for calculating the 
2HDM contribution to the oblique parameters [§]: S, T, U , V, W, X compared to the SM with 
Higgs mass rrtff, as well as the contribution to the muon anomalous magnetic moment. In order 
to take into account collider and flavour limits the program contains interfaces to HiggsBounds 



(version 2) [10] and Superlso [11] respectively. 



7. Usage/parameter-setting/input-output 

The 2HDMC program has been written in a modular fashion with an object oriented structure 
using C++. It can either be used in library mode where the different methods are called by the user 
program or as "ready to compile" command line type programs which are distributed together with 
the main code. 

The code, manual, and full class documentation can be downloaded from http: / /www. 
isv.uu.se/thep/MC/2HDMC. In order to compile the program one needs a gcc compiler 
(3.4 and 4 tested) and the GNU Scientific Library (GSL) installed. There is a template makefile 
which should be adapted to the local conditions and there is a test program Demo with which the 
successful installation can be tested. 

The code consist of four main classes 

• SM: to give the Standard Model input parameters with running a s and quark masses 

• THDM: to specify a general two-Higgs doublet model in terms of the Higgs and Yukawa 
sectors 

• Constraints: to calculate theoretical and experimental constraints on the 2HDM 

• DecayTable: to calculate the decay modes of 2HDM Higgs bosons and the top quark 
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For a more complete description we refer to the class documentation which can be found on the 
2HDMC homepage. 

In order to simplify the communication with other programs there is also a LesHouches style 
input/output format and one can also export a given 2HDM as a user model for MadGraph/MadEvent [ |l2| ] 
to calculate cross-sections etc. 



8. Summary and conclusions 

The 2HDMC program is a versatile and flexible program for exploring the physics of general 
CP-conserving 2HDMs: all the way from defining the model and checking it against theoretical 
and experimental constraints, to exploring its phenomenology such as calculating decay widths 
and scattering cross-sections. Some examples of its use can be found in [|13|, [14], [15, 16 1. 
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